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Abstract 

The magnetic sector of SU(N) Yang-Mills theory at finite temperature is studied. At 
low temperatures, T < 2T C , the analytic expressions for the temperature dependence of the 
magnetic correlator, of the magnetic gluon condensate and of the spatial string tension are 
obtained. Fair agreement with lattice calculations for spatial string tension is obtained for 
SU(2) and SU(3) gauge theories. The relative contribution given by non-zero Matsubara 
modes to the spatial string tension is calculated. At T = 2T C this contribution is of the 
order of 5%. The behavior of magnetic correlator at high temperatures is investigated 
and it is shown that gluo-magnetic condensate increases with temperature as (H 2 )(T) = 
43 ' const • g 8 (T)T 4 . 
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1 Introduction 



It is known that at finite temperature non-Abelian SU(N) gauge theories undergo a de- 
confining phase transition passing from the low-energy glueball phase to the phase of hot 
gluon matter. At the critical point T c where the phase transition occurs the behavior 
of the thermodynamic properties of the system, such as energy density e, specific heat, 
non-ideality (e — 3P)/T 4 , etc., is drastically changed [U 121 El- More than that, the phase 
transition in non-Abelian gauge theories is characterized by the radical rearrangement of 
the non-perturbative (NP) gluon vacuum. 

Recently the behavior of the gauge-invariant two-point correlation functions of the 
gauge fields across the deconfinement phase transition was investigated by numerical 
simulations on a lattice both for the pure-gauge SU(3) theory and for the full QCD 
with two flavors [3J. These data clearly demonstrate strong suppression of the electric 
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component of the correlator above T c and persistence of the magnetic component. The 
magnetic correlator and its contribution to the gluon condensate is kept intact across the 
phase transition temperature, while the confining electric part abruptly disappears above 
T c , so that the electric gluon condensate drops to zero at the deconfining phase transition 
point jlj. These results are completely in line with theoretical predictions of deconfining 
phase transition within the "evaporation model" jSJ approach. Within the framework 
of the effective dilaton Lagrangian at finite temperature the discontinuity of the gluon 
condensate at T = T c was studied in jB] (see also [7j). 

At T = T c the physical string tension becomes zero. It was understood however that in 
non-Abelian gauge theories the Wilson loop for large space-like contours obeys the area 
law at arbitrary temperature. This phenomenon is known as "magnetic confinement" 
and yields non-zero spatial string tension a s {T) jHj. The temperature dependence of 
c s {T) was studied on a lattice for the pure-gauge SU(2) P and SU(3) [21 ED] theories. It 
was shown there that cr s (T) smoothly passes through the phase transition temperature 
and increases with temperature. For temperatures larger than 2T C the scaling behavior, 
^/a s (T)/g 2 T ~ const, settles on [aHUEI]. This regime is imposed by the non-perturbative 
magnetic scale ~ 5 ,2 7 i pi]. It is also firmly established that the 4d SU(N) gauge theory 
at finite temperature is well described by an effective 3d gauge theory with an adjoint 
Higgs field (the so-called dimensional reduction) ^21 US]- Moreover, the scalar Higgs 
fields only weakly influence the physical properties of the gauge theory magnetic sector 
EH ECU UH EE3 CHI- Thus the behavior of certain gluo- magnetic quantities in 3d are 
of direct relevance for the behavior of these quantities in the 4d gauge theory at high 
temperature. In particular, the spatial string tension a s {T) at high T coincides with good 
accuracy with lattice results for the 3d string tension ct 3 jTHl Ej . 

In the present paper we study the magnetic sector of the 4d SU(N) Yang-Mills theory 
at finite temperature. It is demonstrated that thermal behavior of magnetic properties of 
the system is qualitatively different in two temperature regions. We will define T < 2T C 
as low temperature region, and T > 2T C as high temperature region 1 . The temperature 
behavior of the gluo-magnetic correlator, of the magnetic condensate and of the spatial 
string tension in low temperature region (T < 2T C ) are obtained analytically. The relative 
contribution of non-zero Matsubara modes to u s {T) is calculated and found to be of order 
of 5% for pure SU (3) gauge theory at T = 2T C . The temperature interval corresponding 
to the onset of the 3d regime of the Yang-Mills theory is investigated. We consider the 
thermal character of the gluo-magnetic correlator in the scaling region for the spatial string 
tension, ^fo~ s ~ g 2 T. It is shown that the non-perturbative gluo-magnetic condensate 
grows with temperature as ~ g 8 (T)T A at high T. 

2 Gluo-magnetic correlator and spatial string tension 

Gauge- invariant non-local field correlators are of utmost importance for our understanding 
of the NP QCD dynamics (see a recent review [IE] and references therein). The gauge- 
invariant two-point correlators in the Yang- Mills vacuum are defined as EOl IS] 

D^axjx) = (g hrG^jx^jx, 0)G ctA (0)$(0, x)), (1) 

1 "2" should not be considered as an exact number. This matter is discussed in the last section. 
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where = t a G a ^ v is the field-strength tensor, t a are the SU(N) gauge group generators 
in the fundamental representation, tit a t b = \S ab , and 

$(x,y) = Pexp{ig f dx„A„(x)}, A» = t a A^ (2) 

J X 

is the Schwinger parallel transporter with the integration from x to y along a straight- 
line path. These phase factors $ are introduced to take care of the gauge-invariance 
of the correlators. In general form bilocal correlators ("""J are expressed in terms of two 
independent invariant functions of x 2 ,D(x 2 ) and D\(x 2 ) [TTJ | [*"""" [ l2*T]. 

At finite temperature the 0(4) space-time symmetry is broken down to the spatial 
0(3) symmetry and bilocal correlators are described by independent electric and magnetic 
correlation functions. On general grounds one can write for the magnetic correlators the 
following expression 

Df k {x) = ((7 2 tr^(x)$(x,0)^(0)$(0,x)) 



= 5 ik (D H + D«) + (S ik - ^)x 2 ^, (3) 

x z ox z 

where Hi = \eij k Gj k are the magnetic field operators. In a well known way [20] one can 
obtain the area law for large spatial WMoops of size L ^> £ m , where £ m is the magnetic 
correlation length 

(W {C)) sp atiai ~ exp{-a s S min }, (4) 
where the spatial string tension <j s is 

<?s = \J d 2 x(g 2 trH n {x)$(x, 0)# n (0)$(0, x)) (5) 

and H n is the component orthogonal to the space-like plane of the contour. 

Substituting (jBJ) into (|HJ) and keeping in mind that the terms with form the com- 
plete derivative, one can present a s in the following form [2*2*] 



d 2 xD H {x). (6) 



Note that only the term D H enters in the expression ("""} for a s . At T = due to the 0(4) 
invariance electric and magnetic correlators coincide and cr s = a, where a = (420MeV) 2 
is a physical (temporal) string tension. 

Note that magnetic correlator Df k can be expressed through matrix in the adjoint 
representation. Since integration in (J2J) goes along straight-line path, connecting points x 
and y, one has 

D? k (x,y) = (g 2 H?(x)U^(x,y)H b (y)), 
U%{x,y) =2tr (t a $(x,y)t b $(y,x)) . 

The functions D H and contain both perturbative (oc 1/x 4 ) and nonperturbative 
(oc exp{ — |a?|/£ m }) terms and only non-perturbative parts of the correlators contribute to 
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the string tension (for the detailed discussion of these questions see the review paper [TH]). 
The magnetic gluon condensate can be also determined through the NP contributions to 
the correlator and expressed in terms of the functions D and D± at i 2 = 0. The 
gluo-magnetic condensate can be written as 

D*(0) = (ghrH 2 ) = l -(H 2 ) = 3(D H (0) + D*{0)), (8) 

where (H 2 ) = ((gH*) 2 ). 

The lattice data jlj demonstrate that at least in the region of "measured" tempera- 
tures, T < 1.5T C , the NP functions D H = Bq exp{— |a?|/£ TO } and D 1 ^ = B\ exp{ — |a;|/£ m } 
are almost independent of temperature. Correlation length does not change and equals 
to it's value at T = 0, and B\ ps 0.05£> B . Only B slightly grows with T, which 
leads to the slow growth of condensate (H 2 ) at low temperatures. 

Let us see how one can obtain the analytical expressions for the magnetic correlator and 
correspondingly for the spatial string tension and magnetic gluon condensate as functions 
of the temperature. It is well known that the introduction of the temperature for the 
quantum field system in thermodynamic equilibrium is equivalent to compactification 
along the euclidean "time" component £4 with the radius (3 = 1/T and imposing the 
periodic boundary conditions (PBC) for boson fields (anti-periodic for fermion ones). 

Thermal vacuum averages are defined in a standard way 

(...>/»=—/ [DA]...e- s »W (9) 

JPBC 

where partition function is 

Z p = f [DA}e- s ^ A \ S p = f dx 4 [ d 3 xL YM . (10) 

JPBC JO J 

As mentioned above, in the low temperature region gluelump mass M (inverse mag- 
netic correlation length l/£ m ) does not depend on T. Thus we will use zero-temperature 
expression for Df[(x, X4) with PBC to obtain thermal correlator T)f k = (g 2 H^U^H^.) p 
built from gauge fields Hi(x, x±) and phase factor ?7 a dj [A(x, X4)] at a particular "time" 
coordinate £4 

+00 

2>£(x,/?)= D ^,n(3). (11) 

n=— 00 

Usually, this approximation is valid in the case of free noninteracting fields. For 
example, thermal Green function of free scalar field is defined as sum over Matsubara 
modes of propagator at T = 0. At the same time, interaction changes particle mass 
and makes it dependent on T. In our case in the region of low T we take into account 
only "kinematic" change of correlator, and consider it as zero-temperature correlator with 
PBC. Physically the expression (|TT| is valid at /3 > 2£ m , i.e. when the two neighboring 
exponential functions do not overlap. We use the value l/£ m = 1.5 GeV and it is thus 
required that T < l/2£ m = 750 MeV. In the temperature region we consider in the 
present section, T < 2T C ~ 600 MeV, this condition is fulfilled. Moreover comparison of 



4 



calculated below physical quantities with lattice results confirm validity of expression (fTTj) 
at T < 2T C . 

It should be noted that using Poisson's summation formula 

+00 +00 

e ^ x = S ( x ~ *P)> = 27in/(3 (12) 







n=— 00 



one finds, that (fTTJ) is a representation of magnetic correlator D as a Fourier sum over 
Matsubara frequencies 

^A^(x,^) = i^A^(x,^), (13) 



where D(u) = J dr exp(iuT)D(T) is a Fourier image of function D. 
Thus consider the function 

+00 

/(x,/?)= e- M V* 2 +^\ M = l/U- (14) 

n=— 00 

In order to carry out the Matsubara summation over the frequencies uo n = 27rnT in (114}) 
use can be made of the integral representation of the exponent in (fT4"|) . This yields 



/( X> 0) = _1 r^L e -s-M^/,s y e -MW/ 4s _ (15) 

vWo „frt, 

The following summation equation holds 

+00 

E b 2 n 2 Z-V/'/lS \ l:v -.t~ . 

e-TT" = £ e ~ (16) 

n=— 00 

It enables to write (J15|) in the form 



71= — CO 



+ OO „ , +OO „ „ 

b 2 n 2 4vr 2 n 2 



b 

n=— 00 



^ X '^ = M0 E / ^ e " aS "" /S ' ( 17 ) 

n=— 00 ^ 

where a = 1 + 4vr 2 n 2 /M 2 f3 2 and 6 = M 2 x 2 /4. 

The integral in the right-hand side of (|17|) is expressed in terms of the Macdonald 
function K\. Finally we get 

+00 

/(x,T) = 2MT|x| im i Kid^WM' + ul). (18) 

From the asymptotic behavior of function / at short distances we can obtain the tem- 
perature dependence of the magnetic gluon condensate, (H 2 )(T) = (if 2 )/(|x| — > 0,T). 
Using the behavior K\(x — > 0) — > 1/x one has 

\ / rn — — it \ / 
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and at T < M 

(H 2 )(T)/(H 2 ) = 1 + 2e- M ' T + 0{e- 2M ' T ). (20) 

This is completely in line with the lattice result jl] showing the slow increase of the 
magnetic gluon condensate at low T < 1.5.T C . 

Using the expression (jHJ) for a s and equation (fTHj) we get 

! (21) 

n=—oo n 

In arriving to ()21|) use was made of the normalization condition, a s (0) = a, and 
J °° x 2 dxK\{cx) = 2/c 3 . Performing summation in ()21|) we arrive to the expression for the 
temperature dependence of the spatial string tension at low temperatures 

<T) = sinh(M/T) + M/T 

a cosh(M/T)-l ' 1 ' 

and a s (T < M)/a = 1 + 2(M/T + \)e' M l T + 0{e~ 2M ' T ). 

From (|2*T|) we can determine the relative contribution, A n (T) = cr"^°(T)/cr s (T), of the 
non-zero Matsubara modes to the spatial string tension as a function of T. Sorting out 
from (|2*Tj) the term with n = we get 

A (T) - 1 - 4T °osh(M/T)-l (23) 
nl j M sinh(M/T) + M/T 1 j 

It is clear that «7j=°(T) + a^°(T) = a s (T) and A n (0) = 1. 

Also one can obtain from (J2~T|) contribution of individual n-th non-zero Matsubara 
mode 



<(T) M 4 



a 



«=°(T) ^ ( 1 + M 2 /w 2)2 



(24) 



Thus, it can be seen that at T > M/2n contribution of n-th non-zero mode to the spatial 
string tension as compared to the contribution of zero mode, n — 0, very fast "dies out" 
with increasing T 



\2nTj \2ttT ) n 6 + " " 1 &J 



ar°(T) 

3 High temperatures 

According to its definition the quantity M = l/^ m is the inverse magnetic correlation 
length. At T = the correlation lengths were determined by lattice calculations [2*3] . 
On the other hand M can be identified as the mass of the lowest magnetic gluelump 
with quantum numbers J PC = V . Gluelumps are not physical objects and their 
spectrum cannot be measured in experiment. However, they play fundamental role in 
the nonperturbative QCD since gluelump masses define the field correlators in the QCD 
vacuum and in particular the string tension at T = 0. The masses of gluelumps were 
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calculated analytically within the framework of the QCD sum rules QCD string 
model j2H| and computed numerically on the lattice [2Z1- From all the calculations listed 
above one can deduce the value of the correlation length in the interval £ ~ 0.1 4- 0.2 fm. 

In what follows we shall use the value M = 1.5 GeV for inverse magnetic correlation 
length at T = 0. In Fig. 1 we present the dependence of the function A n (T) on T/T c 
(T c = 270 MeV for pure-gauge SU(3) theory). It follows that A„(2T C ) ~ 0.05 , i.e. the 
contribution non-zero (non-static) Matsubara modes into a s is about 5%. This is of the 
same order as the error of the lattice calculations of the spatial string tension [21 El EI] ■ 
Thus at temperatures under consideration, T > 2T C , the main contribution into the 
dynamics of the gauge theory magnetic sector comes from zero (static) Matsubara mode. 
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Figure 1: Relative contribution, a^°/a s , of the non-zero Matsubara modes to the spatial 
string tension for pure-gauge SU(3) theory. 

Next we note that at any T the SU(N) gauge theory is in the phase of magnetic 
confinement. Therefore at high T the exponential form of the magnetic correlator remains 
unchanged. On the other hand the two-point correlation function is determined by the 
amplitude and correlation length which are temperature dependent. Hence one can write 
the magnetic correlator at high temperature T as 



(26) 



The term ~ 0(...), as explained above, contributes neither into gluo-magnetic condensate, 
nor into spatial string tension. Then from © and (|2l)j) we get 



7T 



(27) 



As it was shown in (21 El EI] the temperature dependence of the spatial string tension 
u s (T) at T > 2T C with good accuracy coincides with the behavior of the string tension 03 
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in 3d Yang-Mills theory. The three-dimensional Yang-Mills theory is a superrenormaliz- 
able theory and all physical quantities are determined by the only dimensionful coupling 
constant g 2 = g 2 T. It was shown in [21 El OS] that 



yfijTj = c a g 2 (T)T, (28) 
where use was made of the two-loop expression for g 2 {T) 

9 -(T) = 2M»£ + ^(2>»£), 

\ v a ( 29 ) 
11N , 34 / N 

h 



u 48vr 2 ' 3 \lQn 2 

The two constants c a and A CT were determined using a two-parameter fit to lattice results. 
For the SU(2) gauge theory c a = 0.369 ± 0.014, A CT = (0.076 ± 0.013)T C [9J, while for 
SU(3) gauge theory c a = 0.566 ± 0.013, A CT = (0.104 ± 0.009)T C [2j. In Figs. 2 and 3 we 
present the plots a s (T)/a for SU(2) (Fig. 2) and SU(3) (Fig. 3) gauge theories as functions 
of T/T c [Tc U{2) = 290 MeV, T c SC/(3) = 270 MeV). Solid lines correspond to Eq.® and 
the dashed ones to Eqs. (j28l29|) . The lattice data are from Refe.j2El- 

Further, from (|2*7jl it is possible to find (H 2 )(T). The inverse magnetic correlation 
length l/£ m at high temperature behaves as 

l/U(T) = c m g 2 (T)T, (30) 

where c m is some constant. Using (J2ZJ), flZHl) and (jHOj) we obtain the following expression 
for the temperature dependence of the gluo-magnetic condensate 

(H 2 )(T)=c H g 8 (T)T\ (31) 

and 

c H = -clc 2 m . (32) 

7T 

From (|3^|) one can estimate ch- For instance in SU{2) Yang-Mills theory, c CT ~ 0.37[9 , 
c m ~ 0.92 [2HI, and one finds ch ~ 0.22. However, the value of ch should be found from 
numerical simulations on a lattice for thermal magnetic correlator. 

The result ()31|) seems quite natural from the 3d gauge theory standpoint. As was 
already indicated, the thermal behavior of the physical quantities in the magnetic sector of 
4d SU(N) gauge theory at high T coincides with the behavior of corresponding quantities 
in the 3d gauge theory. Then the NP gluo-magnetic condensate is determined by the 
value of the dimensionful coupling constant #3 and (H 2 )(T) = const g\. 2 Therefore at 
high temperatures (in the scaling region for the spatial string tension) the gluo-magnetic 
correlator has the form 

VgfrT) = l -5 lk c H g\T)T'e-^ T ^. (33) 
6 

The amplitude of this correlator rises as T 4 while the correlation length drops with growing 
temperature as 1/T in the high temperature region. 

2 Strictly speaking, in superrenormalizable 3d Yang-Mills theory vacuum averages either are equal to 
zero, or are defined through corresponding power of coupling constant (73. Disappearance of magnetic 
condensate at high temperatures means that a s = 0, and therefore contradicts lattice data on magnetic 
confinement. 
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Figure 2: Spatial string tension a s (T) / a for SU(2) gauge theory as function of T/T c . Solid 
lines correspond to Eq.(j22J) an d the dashed ones to Eqs. (j28l29|) . The lattice data are from 
Refs-P- 
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Figure 3: Spatial string tension a s (T)/a for SU(3) gauge theory as function of T/T c . Solid 
lines correspond to Eq.(J22J) and the dashed ones to Eqs. (|28|29jl . The lattice data are from 
Refs.|2]. 

4 Discussions and conclusion 

The results which were obtained above allow to estimate the temperature region within 
which the change of the regime in the magnetic sector of SU(N) gauge theory occurs and 
it becomes possible to describe the dynamics in terms of the static correlation functions, 
i.e. in terms of the 3d gauge theory. On physical grounds it is clear that when the 
compactification radius (3 (the inverse temperature 1/T) becomes much smaller than the 
typical dimension I of the system (which in turn is determined either by magnetic string 
thickness, or by the radius of the magnetic gluelump) I = 1/ y/a s (T), the correlators 
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"cease to feel the time coordinate £4" and become pure static. 

Thus the condition T = a/ c s (T) enables to estimate the temperature at which the 
change of the regimes occurs. Consider for defmiteness the SU(3) gauge theory. If one 
starts from the low temperature region then from the numerical solution of the equation 
T_ = a/o" s (T_) and from Eq. (22) for a s {T) one finds T_ = 1.42T C . Alternatively if one 
starts from the high temperature region then the solution of the equation T + = a/ cr s (T + ) 
with a s {T) defined by (|28I29|) gives T + = 2.74T C . Thus we may conclude that within the 
interval T_ < T < T + the change of the regimes takes place and transition to the reduced 
3d Yang-Mills theory occurs. Therefore the commonly accepted value 2T C seems quite 
natural from the outlined physical point of view. 

The growth of the magnetic condensate, (H 2 )(T) oc T 4 , with the increase of the tem- 
perature is thermodynamically advantageous for the quantum field system. The vacuum 
energy density is connected via the trace anomaly 3 with the magnetic condensate by the 
relation e vac = (9 00 ) = — (&/647r 2 )(ff 2 ), b = lliV/3, and the increase of (H 2 ) lowers the 
vacuum energy density. 

In the present paper we have investigated the magnetic sector of SU(N) Yang-Mills 
theory at finite temperature. At low temperatures, T < 2T C , the analytic expressions 
were obtained for the temperature dependence of the magnetic correlator, magnetic gluon 
condensate and the spatial string tension. Fair agreement with lattice calculations for 
cr s (T) was obtained for SU(2) and SU(3) gauge theories. It is demonstrated, that the 
contribution of n-th non-zero Matsubara mode to the spatial string tension as compared 
to the contribution of zero mode, n = 0, very fast "dies out" with increasing temperature 
~ (M/27mT) 4 . The relative contribution given by non-zero Matsubara frequencies to 
the spatial string tension was calculated. At T = 2T C this contribution is about 5%. 
The estimate was presented of the temperature region corresponding to the transition 
to the description of the nonperturbative dynamics in the magnetic sector in terms of 
the 3d correlation functions. We have found the behavior of magnetic correlator at high 
T > 2T C and have shown that gluo-magnetic condensate increases with temperature as 
(H 2 )(T) = c H g\T)T\ 
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